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Abstract 

The transient response of a 0.4 by 0.6 cm rectan- 
gular bore rail accelerator was analyzed using a 
three-dimensional finite-element code. Results are 
presented for the case of a 210-kA current input and 
5-cm arc length. Typically, the copper rail deflected 
to a peak value of 0.08 mm in compression and then 
oscillated at an amplitude of 0.02 mm. Simultaneously 
the insulating side wall of glass fabric base, epoxy 
resin laminate (G-10) was compressed to a peak value 
of 0.13 mm and rebounded to a steady state in 
extension. Projectile pinch or blow-by due to the 
rail {or side-wall) extension or compression, respec- 
tively, can be identified by examining the time history 
of the rail (or side-wall) displacement. For the case 
presented, the effect of blow-by was most significant 
at the side wall characterized by mm-size displacement 
in compression. 

Dynamic stress calculations indicate that the 
G-lO supporting material behind the rail is subjected 
to over 21 MPa - at which the G-10 could fail if the 
laminate was not carefully oriented. Results for a 
polycarbonate resin (Lexan) side wall show much larger 
displacements and stresses than for G-10. Therefore 
the tradeoff between the transparency of Lexan and the 
mechanical strength of G-10 for side-wall material is 
obvious. 

Displacement calculations from the modal method 
are smaller than the results from the direct integra- 
tion method by almost an order of magnitude, because 
the high-frequency effect is neglected. However, this 
effect is significant and dominating for a highly- 
impulsive, wave-propagation problem such as the rail 
accelerator structure. 

Introduction 

A rail accelerator imparts high velocity to a 
projectile through the application of large impulsive 
forces of short duration. For the case of a plasma 
armature, the rail accelerator structure is dynamically 
stressed by the plasma pressure and magnetic forces 
immediately behind the projectile. [IJ The impulsive 
loading causes stress waves to propagate through the 
structure inducing responses in the materials at the 
wave speed. Displacements of the rails and insulating 
side walls ahead of the projectile may hinder the per- 
formance of the accelerator [2] due to a combination 
of two effects: pinch and blow-by. The former occurs 

when the rails or side walls pinch into the projectile 
and frictional interaction intensifies; the latter 
occurs when the clearance between the projectile and 
the rails or side walls opens allowing the plasma to 
blow by the projectile. These effects can be detri- 
mental if the design of the rail and supporting 
structures are insufficient to handle the intensive 
pulse loading of hundreds of megapascals. 

Numerical analysis can be used to estimate the 
structural response of the rail accelerators. For 
some rail geometries, one or two-dimensional analyses 
[3,4] have been applied to estimate the maximum rail 
deflections and the relative effects of different 
structural or material properties. However, it is 
desirable to use a three-dimensional transient program 
to model a rail accelerator structure to visualize the 


transient deflections of the rail in three dimensions 
as well as the propagation of the stress wave. 

The primary purpose of this paper is to predict 
the transient rail and side-wall displacements of a 
0.4 cm by 0.6 cm rectangular bore, one meter long rail 
accelerator using the three-dimensional finite-element 
code MARC on CRAY/IBM computers. Both rail and side 
wall deflections are estimated by applying a pulsed 
plasma loading on both the side wall and the rail and 
a magnetic loading on the rail behind the projectile. 

A secondary purpose is to calculate the stress of the 
supporting materials, including a glass fabric base, 
epoxy resin laminate (G-10), a polycarbonate resin 
(Lexan), and a phenolic resin. 

Results are presented using the Newmark direct 
integration method for a typical case of 210-kA 
current input, with a 5-cm plasma arc length at an 
instantaneous velocity of 2.5 km/s. A modal method 
using 10 modes is also used for comparison. 

Analytical Procedures 

When a material is subjected to a significant 
dynamic loading, disturbances are propagated through 
the body as stress waves. Deformations resulting from 
the impulsive loading will be highly localized within 
the region influenced by the stress wave, while the 
global structure effect is small. For this case wide- 
spectrum frequencies are excited; however, the high- 
frequency modes usually dominate throughout the time of 
interest. [5] 

Finite-element Code 

The three-dimensional finite-element computer 
code, MARC [6], is available in both modal method and 
direct integrating method. The former allows a choice 
of 10 modes and is generally good for low-frequency 
structural dynamic problems. The latter is appropriate 
for high-loading, wave-propagation problems such as 
the rail accelerator. The direct Newmark integration 
method was selected for the analysis, and the 10-mode 
modal method was used for comparison purposes. 

The finite element mesh plot was generated and 
shown in Fig. 1(a). Only a quadrant of the rail ac- 
celerator, of Figs. 1(b) and (c), is needed because of 
two planes of symmetry. The element used is a 8-node, 
isoparametric, three-dimensional hexahedron block. 

The mesh model consists of 300 elements and 504 nodes. 

A typical run on the CRAY-IS uses about 18 minutes (in 
CPU). All the outputs and post plots are then trans- 
ferred to a front computer the IBM 3033. 

For the cases studied, the materials used for the 
rail accelerator were held within the elastic limits. 
Any plastic distortion or non-linear effect is unde- 
sirable. Furthermore, the secondary effect of damping 
is not included in this report which represent the 
upper bound values of both displacement and stress. 

Input Conditions 

A rail accelerator one meter long with a bore of 
0.4 by 0.6 cm was simulated. Rails are made of copper, 
and side-wall and supporting structures are made of 
G-lO/Lexan and phenolic resin. Their properties are 
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listed in Table I. A high current plasma armature was 
produced to accelerate the projectile between the 
rails. 


Referring to Fig. 1, the projectile is set to a 
starting position (reference nodal point 2) and is 
accelerated by the arc plasma covering bore walls of 
element numbers 1 and 241 along the x-direction between 
the rail, all the way through to the muzzle end (refer- 
ence nodal point 21). 

Loading Conditions 

Loading inputs include both arc plasma pressure 
(Pp) and rail magnetic pressure (P^). Pn is 
acting on the back of the projectile of an area (Ap) 
and is propotional to the Lorentz force (Fp). Fp 
is obtained from the input current (I) and^the induct- 
ance gradient (L‘) of the rail: 


L'r 


P 

"^p ■ A„ " 2 h„ S 
P P 


where hp is the rail height, and s is the distance 
between the rails. P^ is equal to the magnetic 
force (Fp) per unit rail area (Ap) and is acting 
on the rail only: 


2 n h 
P 


dL* 

d(S/h) 


where h is the arc plasma height, and dL'/d(S/h) is the 
inductance gradient calculated for this geometry. [7J 

It is assumed that Pp is uniformly distributed 
behind the projectile within an arc length (L^). The 
acting period of time (at) of Pp on a specific ele- 
ment wall is then L^/Vfy^; where is the instan- 
taneous speed of the arc traveling through the bore, 
and La can be estimated by combining available 
measured values [8,9] and the scaling laws. [10] 


Typical values used for the calculations are 
1 isted in Table 1 1 . 


For illustration, the transient pressure inputs 
are plotted in Fig. 2 for rail elements 1 through 4, 
and also for side-wall elements 241 through 244. Where 
a flat-topped pulse of 20 ysec is assumed for Pp of 
450 MPa and a constant Pp is 240 MPa, With known 
arc length and speed, only the arc shape is assumed 
for the transient loading inputs. 

Different pulse shapes of inertia loading and re- 
sponses have been well investigated. [11] The square 
wave renders the maximum response, by an amplification 
factor of two, as compared to other shapes: sinus- 

oidal, triangular, sudden-rise, or sudden-fall. The 
flat-topped shape pulse assumed in the paper should 
give more conservative values of displacement, although 
the actual shape of loading may be quite complicated. 

Results and Discussions 


The primary case (P41) is calculated using G-10 
as the side-wall insulator and applying the direct 
integration method. Two secondary cases are included 
for the relative comparison of using transparent Lexan 
as the side-wall material (Case N41) and of applying 
the modal method with ten modes (Case Q41). 


G-10 backing material behind the rail, and the sec- 
tional stress contours of case P41. 

Figure 3 shows the rail displacements at the 
breech end (nodes 1 and 2) and the muzzle end (nodes 20 
and 21). Each nodal point has a similar response 
characterized by a pulse due to the impulsive loading 
and then approaching a quasi-steady value. Figure 4 
shows that the typical response of a single node 3, 

5-cm downstream of the projectile's starting position 
which is another 5 cm from the breech end. It also 
shows the peak displacement reaches 0.06 mm of 
compression and oscillates, with a period of 30 usec 
within 0.01 to 0.03 mm range. The instantaneous 
position of the projectile corresponds to a 40 psec 
time scale in Fig. 4, where a compression of 0.03 mm 
indicates possible blow-by which means the plasma 
leaks through the clearance to the front of the 
projectile. The actual processes of blow-by or pinch 
are complicated because of their high-speed, nigh- 
temperature, and dynamic features. However the trans- 
ient displacement can be used as a critical parameter 
to assess the possibility of blow-by or pinch. Crite- 
rion for the parameter may be different for different 
accelerators under various operating conditions. 

Corresponding side-wall displacements are plotted 
in Fig. 5, at the breech end (nodes 421 and 422) and 
at the muzzle end (nodes 440 and 441). Peak deflec- 
tions of between 0.07-0.15 mm are obtained along the 
one meter length of the accelerator, while the quasi- 
steady deflections are under extension (indicated by 
minus values). Figure 6 shows the displacements of 
node 423. At the current projectile position 5-cm 
downstream (corresponding to 40 psec), a compression 
of 0.07 mm indicates that the effect of blow-by may be 
larger than that on the rail. However, the pinch 
effect seems to be negligible on the projectile 
because the displacements are mostly positive at that 
1 nstant. 

The transient stress calculation for G-10 element 
83 (which is under element 23 in Fig. 1(a)) exhibits a 
typical variation during the time of interest as shown 
in Fig. 7. This element behind the rail, at 5-cm 
downstream of the projectile's starting position, has 
a stress wave form similar to these of displacement 
response, peaked at 22 MPa and approached quasi-steady 
value of 8 MPa. If the G-10 laminates are oriented in 
the weakest direction at which the yield strength is 
approximately 7 MPa, failure can occur. [12] Stress 
contours are plotted at a section cut through the rail 
and its G-10 backing material (elements 61 through 100 
and 241 through 260), for each 100 psec time period. 

The traveling characteristic of the stress wave can be 
seen from Fig. 8 where the stress concentration reaches 
410 MPa on the rail for most of the time while peaked 
at 700 MPa at the end due to rebounding at the muzzle 
end. 

Secondary Case (N41) 

It can be seen from Fig. 9 that the deflection 
reaches 0.8 mm using Lexan as the side-wall material. 
This displacement is about an order of magnitude larger 
than the G-10 case shown in Fig. 6, and the effect of 
blow-by will be significant. The use of Lexan at high 
current (210-kA) seems to be unfeasible. The trade 
off between the transparency and the operating cur- 
rent, or the mechanical strength, is thus obvious. 

Secondary Case (Q41) 


The Primary Case (P41) 

Time history results are presented for both rail 
and side-wall displacements, the element stress of 


For the purpose of demonstrating the difference 
between the application of numerical methods, the same 
input conditions of the primary case is calculated by 
using the modal method. As shown in Fig. 10, rail 
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displacements are much smaller. On the other hand 
periods are larger when compared with Fig. 4. The 
reason for these differences is the modal method does 
not include enough modes so that the dominating effects 
of high frequencies of the wave are neglected. On the 
other hand, the direct integration method is appro- 
priate for the highly impulsive, wave-propagation 
problem such as the rail accelerator. 

Concluding Remarks 

The response of a rail accelerator structure to 
an impulsive loading can be predicted by a three- 
dimensional, finite-element code - MARC using the 
direct integration method. The effects of blow-by or 
pinch are shown by rail or side-wall displacements. 
Mechanical limits are also discussed to the first order 
by calculating the dynamic stress values. Relative 
merits of different materials or geometry can readily 
be assessed. 
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TABLE I. - PROPERTIES OF RAIL ACCELERATOR MATERIALS 



Young's 

modulus, 

GPa 

Poisson 

ratio 

Density, 

gm/cm^ 

Yield 

point, 

MPa 

Copper 

117 

0.37 

8.4 

400 

G-10 Laminate 

17.2 

0.33 

1.75 

7-350* 

Lexan 

2.3 

0.33 

1.2 

63 

Phenolic Resin 

8.3 

0.33 

1.4 

45 


*Depends on orientation of the laminates. 


TABLE II. - VALUES USED FOR CALCULATIONS 


Current 

Arc 

Velocity 

Inductance 

Inductance 

Rail 

Arc 

Rai 1 

I. 

length. 

Vm. 

gradient 

L', 

gradient 

height 

height 

width 

KA 

•-a* 

km/s 

dL-/d(S/h). 

h. 

"p 

s, 


cm 


ph/m 

pn/m 

cm 

cm 

cm 

210 

5 

2.5 

0.48 

0.245 

0.63 

0.38 

0.63 







































(a) Mesh plot of the 
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SYMMETRIC A QUADRANT BETWEEN 
PLANE B TWO SYMMETRIC PLANES 

I IS USED FOR MODELING 



(b) Rail accelerator cross section. 


Figure 1. - Continued. 



(c)Rail accelerator (gun). 
Figure 1. - Concluded. 
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Figure 2. - Transient profiles of pressure loading for the 
first four elements of both the rail (1-4) and side wall 
(241-244). 
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Figure 3. - Rail displacements at breech and muzzle ends. 
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Figure 4. - Rail displacement 5-cm downstream. 
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Figure 5. - Side-waii displacements at breech/muzzle ends. 
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Figure 7. - Stress in G-10 backing material 5-cm downstream. 
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Figure 8. - Transient stress contours for copper rail (elements 61-80), G-10 backing (elements 81-100), and 
G-10 siding (elements 241-260). 
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